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Abstract 

Objectives: Our aims were to replicate some previously reported associations of single nucleotide polymorphisms (SNPs) in 
five genes (A2BP1, COGS, GDF5, HFE, ESRl) with hand osteoarthritis (OA), and to examine whether genes (BCAP29, DI02, 
DUS4L, DVWA, HLA, PTGS2, PARD3B, TGFBl and TRIBl) associated with OA at other joint sites were associated with hand OA 
among Finnish women. 

Design:\Ne examined the bilateral hand radiographs of 542 occupationally active Finnish female dentists and teachers aged 
45 to 63 and classified them according to the presence of OA by using reference images. Data regarding finger joint pain 
and other risk factors were collected using a questionnaire. We defined two hand OA phenotypes: radiographic OA in at 
least three joints (ROA) and symptomatic DIP OA. The genotypes were determined by PCR-based methods. In statistical 
analysis, we used SNPStats software, the chi-square test and logistic regression. 

Results: Of the SNPs, rs716508 in A2BP1 was associated with ROA (OR = 0.7, 95% CI 0.5-0.9) and rs1800470 in TGFBl with 
symptomatic DIP OA (1.8, 1.2-2.9). We found an interaction between ESRl (rs9340799) and occupation: teachers with the 
minor allele were at an increased risk of symptomatic DIP OA (2.8, 1.3-6.5). We saw no association among the dentists. We 
also found that the carriage of the C0C5 rs3757713 C allele increased the risk of ROA only among women with the BCAP29 
rsl0953541 CC genotype (2.6; 1.1-6.1). There was also a suggestive interaction between the HFE rsl79945 and the ESRl 
rs9340799, and the carriage of the minor allele of either of these SNPs was associated with an increased risk of symptomatic 
DIP OA (2.1, 1.3-2.5). 

Conclusions: Our results support the earlier findings of A2BP1 and TBGFl being OA susceptibility genes and provide 
evidence of a possible gene-gene interaction in the genetic influence on hand OA predisposition. 
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Introduction 

Osteoarthritis (OA) shows clinical heterogeneity in localization 
and progression [1]. Patients may have only one affected joint 
(monoarthritis) at the time of diagnosis, several affected joints 
within a single region {e.g. in hand OA), or several involved joints 
at various sites, e.g. the hip, knee, and hand (polyarticular OA, 
generalized OA). Twin and family studies have demonstrated a 
significant contribution of genetic factors that account for up to 
half of the risk of developing OA [2,3] . 

The hand is among the most prevalent sites affected by OA, 
especially among women over the age of 50 [4]. The simultaneous 
involvement of multiple hand joints makes hand OA a heteroge- 
neous disorder that is complex to study [5]. Hand OA, largely 
mirroring the generalized OA variant, is thought to be more 
heritable than hip or knee OA [3,6]. It is generally accepted that 
as a complex disorder, the development of hand OA is modulated 
by many genes with small effects and gene-environment interac- 



tion. The genetic influence may involve either a structural defect 
(e.g., in collagen), alterations in the structural extracellular matrix 
(ECM) proteins of cartilage and bone, an enhanced inflammatory 
component in the disease process, or a genetic influence on a 
known risk factor for OA, such as obesity [7]. 

Genome-wide scans with different hand OA phenotypes 
suggested that chromosomes 1, 2, 3, 4, 7, 8, 9, 11, 13, 15, 16, 
19, 20 may harbor susceptibility genes [8-13]. Numerous 
candidate gene studies have been carried out to assess the 
association of a particular variant with hand OA. However, 
according to a systematic review, specific associations between a 
gene and hand OA have rarely been analyzed by more than one 
study, and only for two genes [AGCl and HFE) have significant 
associations been replicated by at least two independent studies 
[14]. 

A genome-wide association study (GWAS) is a promising tool 
for discovering the genetic basis of common diseases [15]. GWAS 
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were successful in the identification of 11 loci associated with 
different OA phenotypes, in particular knee and hip OA [16]. 
Nevertheless, under the strict criteria of replication and a known 
functional role, the growth and differentiation factor 5 (GDF5) is 
the only truly OA-associated gene at present [17]. To our 
knowledge only two GWAS of hand OA have been published so 
far [18,19]. 

The first GWAS study detected and replicated an association 
with an SNP (rs7 16508) in the first intron of the ataxin 2-binding 
protein 1 gene (A2BP1) [18]. In the second GWAS, a novel 
common variant (rs3815148) in intron 12 of the component of the 
oligomeric Golgi complex 5 gene (COGS) was associated with 
knee and/ or hand OA in both discovery and replicated samples. 
This SNP is in almost complete hnkage disequilibrium with 
rs3757713 (68 kb upstream of GPR22), which could be the link to 
OA association. In the same GWAS, three loci in the GDF5 gene 
(rs491 1494, rs6088813, rs6087705) were associated with hand OA 
in the discovery sample, but were not replicated [19]. 

In order to devise future prevention and treatment, strategies for 
OA replication studies that verify positive findings from GWAS 
are needed. Such studies will also enable us to determine whether 
the effect is specific to a certain OA phenotype, and/or to a 
particular population. Furthermore, GWAS have made it evident 
that most of the genetic risk of complex diseases will be attributed 
to numerous genes with small to moderate effect sizes [16]. In 
order to understand the manner in which the individual genes that 
are implicated in OA exert their effect, gene-gene combination 
effect and interaction need to be explored. 

Our aims were to replicate some previously reported associa- 
tions with hand OA for single nucleotide polymorphisms (SNPs) in 
five genes (A2BP1, COG5, GDF5, HFE, ESRl) and to examine 
whether variants in nine genes (BCAP29, DI02, DUS4L, DVWA, 
HLA, PTGS2, PARD3B, TGFBl and TRIBl) are associated with 
hand OA among Finnish middle-aged women. All selected SNPs 
were locatc'd on genes with prior evidence from GWAS or 
candidate gene association studies of an association with different 
OA phenotypes. 

Participants and IVIethods 

Study Population 

The study participants were identified from the registers of the 
Finnish Dental Association and the Finnish Teachers Trade 
Union, and randomly selected from both occupational groups by 
using the place of residence (Helsinki or its neighboring cities) as 
an inclusion criterion. The samples were restiicted to women aged 
45 to 63. Of those who received the questionnaires in 2002, 295 
(68%) dentists and 248 (57%) teachers participated in a clinical 
examination. 

Ethics Statement 

Participation in the study was voluntary and based on written 
informed consent. The study proposal was approved by the 
Hospital District of H(;lsinki and Uusimaa Ethics Committee for 
Research in Occupational Health and Safety. 

Hand Radiography and Image Analysis 

Both hands of the study participants were radiographed by 
exposing Kodak X-ray films with Siemens X-ray equipment 
(48 kV, 1 0 mA, focus film distance =115 cm; Siemens, Munich, 
Germany). The analogue radiographs were evaluated by an 
experienced radiologist who was blinded to the occupation, age, 
and all health data of the participants. Each distal interphalangeal 
(DIP), proximal interphalangeal (PIP), thumb interphalangeal, and 



metacarpophalangeal (MCP) joint of both hands was graded 
separately, and classified for the presence of OA using a modified 
KeUgren and Lawrence system [20]. The classification criteria 
were: grade 0 = no OA, grade 1 = doubtful OA, grade 2 = mild 
OA, grade 3 - moderate OA, and grade 4 = severe OA. Reference 
images were used; their description is given elsewhere [21]. A 
second experienced radiologist interpreted 46 randomly chosen 
radiographs. A second reading of these 46 radiographs was 
independendy performed by the primary radiologist (TV). The 
reliability of the OA classification was estimated by measuring 
intra-observer and inter-observer agreements using the weighted 
Cohen's kappa coefficient with quadratic weights [22]. The inter- 
observer agreement for OA ranged from 0.67 to 0.85 (good) for 
DIP joints, from 0.39 to 0.61 (moderate) for PIP joints, and from 
0.18 to 0.69 (poor to good) for MCP joints. The intra-observer 
agreement for OA ranged from good to very good, 0.73 to 0.88 for 
DIP joints, from 0.67 to 0.92 for PIP joints, and from 0.59 to 1.0 
for MCP joints. [21]. 

Questionnaires and Interviews 

All study participants received a self-administered questionnaire 
that included rjuestions on anthropometric measures. Information 
on symptoms (pain, tenderness) in each joint studied was collected 
with the prompt: Tlease point out on the picture below in which 
finger joint you have felt pain or tenderness during the past 30 
days.' The participants were also asked to mark the intensity of the 
pain: 0 = no pain, 1 = mild pain, 2 = moderate pain, 3 = severe 
pain. Weight was measured without shoes to the accuracy of 
0.1 kg. Body mass index (BMI) (weight (kg) per height squared 
(m^)) was calculated on the basis of self-reported height and 
measured weight. 

Hand OA Phenotypes 

Two hand OA phenotypes were used. If the participant had 
radiographic findings (grade &2) in at least three finger joints she 
was classified as having radiographic OA (ROA). Otherwise, the 
participant was classified as not having ROA. If the participant 
had both radiographic findings (grade >2) and symptoms (grade 
&1) in at least two DIP joints, she was classified as having 
symptomatic DIP OA. Otherwise, the subject was classified as not 
having symptomatic DIP OA. 

SNP Selection 

We aimed primarily to replicate the associations for the 
candidate genes that have been identified by two recent GWAS 
of hand OA [18,19] and genes with a known association with hand 
OA patiiology. We chose two SNPs (rs7 16508 and rs3815 148) tiiat 
reached a genome-wide level of significance for association with 
hand OA [18,19], and variants in the GDF5 (rsl43383) [23], HFE 
(rsl 799945) [24], and ESRl (rs2234693, rs9340799) [25,26] 
genes. The latter genes were chosen on the basis of the rapidly 
increasing prevalence of hand OA in women over the age of 45. In 
addition, we searched PubMed for studies reporting an association 
between any OA phenotype and the candidate genes located on 
the chromosomes identified by genome-wide linkage studies as 
harboring susceptibility genes for hand OA [8-13] (Leppavuori 
et al. 1999, Demissie et al. 2002, Stefansson et al. 2003, Hunter 
et al. 2004, Greig et al. 2006, Livshits et al. 2007). Studies 
published by 10.1 1.2010 were reviewed in order to select relevant 
SNPs. Whenever at least one significant functional SNP was 
reported for a given candidate gene in hand OA, the SNP was 
selected for analyses in our samples. Finally, we selected variants 
from nine candidate genes (BCAP29, DI02, DUS4L, DVWA, 
HLA, PTGS2, PARD3B, TGFBl and TRIBl) from seven 
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publications [27-33] for exploratory analysis in our study. Several 
predicted functional effects were found in the F-SNP database for 
the studied polymorphisms [34]. The description of SNPs selected 
for replication is given in Table SI. 

Genotyping Analysis 

Blood samples were taken from each study participant at the 
clinical examination and stored at +4°C until we extracted DNA 
using a DNA extraction kit (PUREGENE DNA Purification Kit; 
Centra Systems, Plymouth, MN, USA). 

The ESRl PvuII, and Xbal polymorphisms (rs2234693, and 
rs9340799 respectively) were genotyped using the RFLP method 
as essentially described in [35]. The TGFBl LeulOPro (29T>C, 
rs 1800470 formerly known as rs 1982073) polymorphism was 
genotyped by the TaqMan PGR method as described in [36] with 
somewhat altered PGR conditions (2 min +50°C, 10 min +95°C 
and 40 cycles 15 s +95°C, 1 min +62°G). 

We analyzed the rest of the studied polymorphisms in A2BP1 
(rs716508), BGAP29 (rsl0953541), COGS (rs3757713 and 
rs3815148), PTGS2 (rs4140564), DI02 (rs225014), DUS4L 
(rs4730250), DWA (rs7639618), GDF5 (rsl43383), HFE 
(rsl799945), HLA (rsl0947262), PARD3B (rsl207421) and TRIBl 
(rs45 12391) simultaneously using OpenArray equipment and the 

TaqMan SNP Genotyping Assays (Applied Biosystems, G 122 

4093_10, G_2618842_20, G_274751 19_10, G_259941 14_10, 
G_31274663_20, G_15819951_10, G_32373604_10, G_11767 
13_10, G_1270479_l, G_1085600_10, G_32201424_10, G_88 
07483_10 and C 310264_20, respectively). 

Genotyping completion rate was from 99.8 to 100%. 

Statistical Analysis 

The potential deviation of the allele frequencies from the 
Hardy-Weinberg equilibrium (HWE) was tested from total 
population using the chi-srjuare test. The degree of pairwise 
linkage disequilibrium (LD) for two SNPs in the COGS gene and 
two SNPs in the ESRl gene were calculated using SNPStats 
software [37]. Each SNP was analyzed in turn. Haplotypes were 
constructed and analyzed by the SNPStats program. Logistic 
regression analysis was used to test the associations between 
selected SNPs and two hand OA phenotypes. For each SNP, a log 
additive model of inheritance was fitted. Gene-occupation 
interaction was tested for all SNPs, to evaluate whether the 
association between the SNPs and OA was modified by 
occupation. In addition, gene-gene interactions and gene-gene 
combination effects were evaluated by a logistic regression model 
with a dummy variable (0, 1) for the SNPs to compare the 
magnitude of their odds ratios (ORs). From each gene, we selected 
a representative SNP that showed the lowest p-value in the 
replication analysis. We used the dominant model, with the 
homozygous genotype of the major allele as the reference. Both 
crude and adjusted ORs and their 95% confidence intervals were 
calculated. We adjusted the ORs for the potential confounding 
effects of age (continuous), occupation (dentists vs. teachers) and 
BMI (continuous). Since the crude and adjusted ORs did not differ 
significantly, only the adjusted ORs are shown. 

All analyses were hypothesis driven. We have used p<0,05and 
p<0.008 as significance level in replication and exploration 
analyses, respectively. In exploratory analyses p-values were 
adjusted for multiple testing using Sidak's method [38]. We used 
SNPStats software [37] for the analyses. 



Results 

The prevalence of ROA and symptomatic DIP OA were 
29.5%, and 9.0%, respectively. The ROA prevalence was 
statistically significantly higher among the teachers than the 
dentists (Table 1). The genotype distributions of the selected SNPs 
did not deviate from the HWE. There was no difference between 
the minor allele frequencies for any SNPs of the two occupations. 

Replication Analyses of Hand OA-associated SNPs 

Of the two SNPs identified in GWAS, only the rs7 16508 located 
in the A2BP1 gene was associated with ROA (OR = 0.68, 95% GI 
0.50-0.93) (Table 2). However it was not associated with 
symptomatic DIP OA (Table 3). The risk of symptomatic DIP 
OA was marginally, though not statistically significantly, elevated 
with the occurrence of the minor alleles of two SNPs (rs3757713, 
rs3815148) in the COGS gene. 

We found no statistically significant associations between the 
SNPs in tiie GDF5, HFE and ESRl genes and hand OA. Linkage 
disequilibrium was very strong between the two COGS SNPs 
(rs3757713 and rs3815148) as weU as die two ESRl SNPs {PvuII 
rs2234693 and Xbal rs9340799) (D' = 0.95, p<0.0001 and 
D' = 0.99, p<0.0001, respectively). Haplotypes comprised of the 
two COGS SNPs or of the two ESRl SNPs were not associated with 
hand OA. 

We found an interaction between the ESRl rs9340799 SNP and 
occupation (p = 0.0064) in relation to symptomatic DIP OA. The 
teachers with the minor allele were at an almost three-fold 
increased risk of .symptomatic DIP OA (OR = 2.84, 95% GI 1.25- 
6.48) while the dentists carrying the allele were at a lower risk 
(OR = 0.50, 95% CI 0.19-1.33) than tiiose with the major allele. 
Adjustment for the use of hormone therapy did not affect the 
observed associations. 

Exploratory Analyses of Nine OA-associated SNPs 

Of nine SNPs associated with different OA phenotypes (knee, 
hip, spine and multiple joints) in previous studies, only the TGFBl 
LeulOPro (rs 1800470) SNP was associated with symptomatic DIP 
OA (Table 3, OR= 1.84, 95% GI 1.16-2.91). We observed a 
suggestive (p = 0.047) interactive effect between the PTGS2 
rs4 140564 SNP and occupation in relation to ROA. The teachers 
with the minor allele were at a 1.44-fold (95% CI 0.69-3.02) 
increased risk, and the dentists were at a lower risk (OR = 0.40, 
95% GI 0.13-1.21) (Table 2). 

Gene-gene Combination Effects and Interaction 

To evaluate possible gene-gene combination effects and 
interactions, we selected two SNPs from the genes located on 
chromosome 7 {COGS rs3757713 and BCAP29 rsl0953541) and 
two SNPs from the genes located on chromosome 6 [HFE 
rsl 79945 and ESRl rs9340799) that showed relatively stronger 
associations with our hand OA phenot^-pes. We found a 
statistically significant interaction (p — 0.034) between the COGS 
rs3757713 and BCAP29 rsl0953541 SNPs and a suggestive 
interaction (p = 0.096) between tiie HFE rsl 79945 and ESRl 
rs9340799 SNPs (Table 4). Carriage of the COGS rs3757713 G 
allele increased the risk of ROA only in women with the BCAP29 
rsl0953541 GG genotype. However, the likelihood of ROA in the 
carriers of the minor allele of either SNPs (rs3757713 or 
rsl0953541) was 1.49 (95% CI 1.07-2.06) times higher than tiiat 
in the non-carriers of these alleles. The carriage of the HFE 
rsl 79945 G allele increased the risk of symptomatic DIP OA in 
women homozygous for the ESRl rs9340799 major allele (AA 
genotype). Carriage of the minor allele of either SNP (rsl79945 or 
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Table 1. Selected characteristics of study population. 







Radiographic OA* 




Symptomatic DIP OA** 




No 


Yes 


No 


Yes 


All (n = 542) 


382 (70.5%) 


1 60 (29.5%) 


493 (91.0%) 


49 (9.0%) 


Age (mean ± SD) 


53.0±5.2 


56.3 ±4.7 


53.6±5.2 


57.7±4.2 


BMI (mean ± SD) 


24.3±3.5 


25.0±3.8 


24.4 ±3. 6 


25.1 ±3.3 


Pinch grip strength(mean ± SD} 


55.9±7.4 


52.3±9.2 


55.2±8.1 


51.5±7.0 


Dentists (n = 294) 


222 (75.5%) 


72 (24.5%) 


274 (93.2%) 


20 (6.8%) 


Age (mean ± SD) 


52.6±5.6 


57.1 ±5.6 


53.3±5.8 


59.4±4.8 


BMI (mean ± SD) 


23.8±3.2 


24.5 ±3.3 


23.9±3.3 


23.9±3.0 


Pinch grip strength(mean ± SD) 


56.3 ±7.5 


52.7 ±10.1 


55.8±8.2 


50.4±8.4 


Teachers (n = 248) 


160 (64.5%) 


88 (35.5%) 


219 (88.3%) 


29 (11.7%) 


Age (mean ± SD) 


53.5±4.5 


55.7±3.8 


53.9±4.5 


56.5±3.3 


BMI (mean ± SD) 


25.0±3.8 


25.4±4.1 


25.0±4.0 


25.9±3.3 


Pinch grip strength(mean ± SD) 


55.3±7.1 


52.0±8.4 


54.4±8.0 


52.2±5.8 



*p = 0.006 for difference between prevalence among dentists and among teachers. 
**p = 0.05 for difference between the prevalence among dentists and among teachers. 
doi:l 0.1 371 /journal.pone.009741 7.t001 



rs9340799) was associated with an approximately two-fold 
increased risk of symptomatic DIP OA (OR = 2.12, 95% CI 
1.28-2.50). 

Discussion 

In this study we attempted to replicate some previously reported 
genetic associations with hand OA. Of the two SNPs identified in 
GWAS, we found a statistically significant association between the 
A2BP1 rs7 16508 SNP and radiographic hand OA. In addition, we 
observed evidence of an interactive effect between the COGS 
rs3757713 and the BCAP29 rsl0953541, and radiographic hand 
OA, and between the HFE rs 179945 and the ESRl rs9340799 
and symptomatic DIP OA. Moreover, our exploratory analysis of 
the SNPs in nine selected genes that have a hypothetical 
association with hand OA provided evidence of an association 
between the TGFBl LeulOPro (rs 1800470) SNP and symptomatic 
DIP OA. 

In order to succeed in the replication of genotype-phenotype 
associations, both the initial and replication studies should 
preferably have the same gene variants, disease phenotypes, 
genetic models, and population segment {e.g. the same gender), and 
should be adjusted for the same confounding factors. We selected 
two hand OA phenotypes (ROA in at least three finger joints and 
symptomatic DIP OA), that are alike or very similar to those in the 
discovery studies. Both our phenotypes were defined on the basis 
of a K— L score of grade 2 or more (at least mUd OA), which has 
been commonly used for classification in other studies [24,26,39]. 
Radiographic hand OA in our study was based solely on 
radiographic evidence of joint destruction, whUe symptomatic 
DIP OA was defined on the basis of the simultaneous presence of 
radiographic changes and pain in the most commonly affected 
joint (DIP). 

GWAS with two discovery and four replication samples 
reported an association between the rs7 16508 SNP located in 
the first intron of the A2BP1 gene and hand OA in the population 
of European ancestry [18]. Meta-analysis suggested that the minor 
allele of this SNP diminished the risk of hand OA by 33-41 %, with 
the protective effect larger for the severe hand OA phenotype (at 
least three joints affected). In line with the initial study, we 



observed a 32% reduction in the risk of radiographic OA per 
minor allele. However, we found no association between the 
rs7 16508 SNP and symptomatic DIP OA. 

The A2BP1 gene is known to be associated with several human 
diseases, such as autism [40], cancer [41] and obesity [42], though 
its role in the development of hand OA needs further exploration. 
Since the A2BP1 gene is abundantly expressed in the skeletal 
muscle [43], it is possible that the association between the 
rs7 16508 SNP and hand OA could be via muscle strength. We 
have earlier reported an association between pinch grip strength 
and symptomatic ROA, but radiographic findings or pain were 
not associated with pinch grip strength [44] . In the present study, 
the rs7 16508 SNP was not associated with pinch grip strength and 
the odds of ROA remained unchanged after adjustment for pinch 
grip strength (data not shown). 

Although we failed, in the sample of Finnish female dentists and 
teachers, to replicate the associations of the other SNPs with 
reported effects on hand OA, the effect sizes of the risk alleles on 
ROA and symptomatic DIP OA were similar to those found in the 
initial studies. Reasons for the lack of replication could be due to 
differences in phenotypes, in minor allele frequencies, and the high 
likelihood of false negative findings, because our study was not 
large enough to detect very small effects. 

Multiple susceptibility variants with small effect sizes may also 
interact with each other or with environmental factors in their 
influence on hand OA predisposition. 

In agreement with this, our result suggests that the effect of the 
rs3757713 SNP (that is in complete linkage with rs3815148) on 
hand ROA might be diluted by the carriage of the BCAP29 
rsl0953541 minor aUele. The C0G5 and BCAP29 genes together 
with the other four genes {PRKAR2B, HBPl, GPR22 and DUS4L) 
are located on chromosome 7q22 within a block with high-linkage 
disequilibrium. The GWAS meta-analyses of knee OA showed 
that this LD block contains an OA susceptibility locus [27]. Both 
genes were expressed in OA joint tissue but were down-regulated 
in OA cartilage compared to normal control cartilage [45]. 

In addition, we found that the effect of the //Fi? rs 17 99945 on 
symptomatic DIP OA was diluted by the carriage of the ESRl 
rs9340799 minor allele. Both genes may play a role in hand OA, 



PLOS ONE I www.plosone.org 



4 



May 2014 | Volume 9 | Issue 5 | e97417 



Replication Study of Osteoarthritis Candidate Genes 



> 
a 



< 
O 



> 
a. 



to 
Q- 



-Q 

c 
to 



< 

O 

cc 



5 

l/l 

— 

Q. 
O 

E 

o 
a. 

-a 

OJ 

t: 



(N 

n 



a 

3 

o 



a 



ca 
u 
z 



01 

V} 



as 


in 


CO 


o\ 


<N 


IT) 






CO 






OJ 


00 




1^ 




rsi 


CO 




CO 




OJ 


00 


m 


1^ 


1^ 


o 


CO 






00 


OJ 




o 


ro 


o 










u-i 




LO 


LO 










CO 




Ln 


o 






O 


on 


CO 




CO 


Ln 








00 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


m 




ro 




LO 






OJ 


rsi 


ro 








00 




m 




Ln 


Ln 










OJ 




m 


rsi 


ro 




OJ 


VO 


in 




p 




00 




LO 




vD 


fN 










00 


OJ 




CO 












04 


p 




rn 


CO 


vq 




in 


Ln 


o 


















1 












1 


1 




















1 






oj 




4- 


ro 


o 








O 


O 


■4- 


CO 


\o 


rn 


4- 






rsi 






■4- 




o 


on 


Ov 


>i> 


a\ 


Ln 






4- 


fN 


av 


m 


ro 


CO 




LO 


00 


CO 


U-i 








1^ 




00 


1^ 




vO 






OS 




LO 








in 






CO 






in 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


















































s 
















































































d 


d 






d 






d 


d 


d 






d 




d 


d 






d 






d 


d 








d 






d 
















OJ 






U-I 






ro 






o 






o 








in 






o 






Ln 






CO 








o 












u-i 






Ln 












rsi 








m 






p 












vO 








d 






d 






d 






d 






d 






d 








d 












d 






d 








CO 


00 


CO 


00 


CO 


CO 


ro 




rsi 


OJ 


rN 


m 






o 


Ln 




ro 


























ro 


ro 
































rsi 


rsi 


rsi 














en 












































































ro 






fN 














^ 














fN 


fN 


fN 




\o 


m 




















LO 






c^ 


CO 






CO 


ro 


ro 


rsi 


ro 


rj 




o 




o 


Ln 














o 


o 


o 




CO 






in 


ro 


ro 


















■=1- 












■ba- 


rsi 


rsi 


rsi 




1^ 


\o 


00 


















sts 


lers 




Sts 


lers 




Sts 


lers 




Sts 


lers 




Sts 


lers 




sts 


lers 




Sts 


lers 






Sts 


lers 




sts 


lers 




Sts 


lers 




Sts 


enti 


u 




enti 


ro 




enti 


u 
ro 




enti 


u 

fD 




enti 


ro 




enti 


u 
ro 




enti 


u 
ro 






enti 


u 
ro 




enti 


u 
ro 

OJ 




enti 


ro 




enti 






< 


Q 




< 


Q 




< 






< 


Q 




< 


Q 




< 








< 


Q 




< 


Q 


1— 


< 






< 


Q 












OJ 












CO 












o 




















00 




















OJ 






ro 


















ro 








00 






CTv 












VO 














a\ 






CTi 


















vD 














m 






lO 
































00 






? 






rsi 




















OJ 














































































ro 






00 












Ln 






m 


























CO 


























ro 

CO 






0^ 






c^ 

v£) 




















04 






vO 






OJ 














IT) 












0^ 






^3" 






o 








m 
o 












0v 














LO 












ro 






0^ 






m 




















O 






ro 






0^ 








rv. 






00 


















rN 






ro 














rsi 






\o 






O 








ro 






rn 












"vi 






OJ 


























VI 


























O 
































































fN 
























« 














ro 


















































w 
































































> 
































































































































C 
































































IB 
































































& 
































































0 
































































to 




























LO 






IT) 






U-1 
























0 


OJ 

LO 






m 
Q 






< 














<J 












Ll_ 






LU 


















a 








cr 












< 








o 






o 






Q 












LO 
LU 






in 

LU 






X 


J5 






< 




















u 






u 






<J 






X 














u 


Q- 












Q 






X 





PLOS ONE I www.plosone.org 



5 



May 2014 | Volume 9 | Issue 5 | e97417 



Replication Study of Osteoarthritis Candidate Genes 



CTi CO CO CO 



CO vO 0^ 



U-i CO CO CO 



< 

O 



I— m 



> 

a 



tN <— -— t— 



tN <— tN t— 



< 



rsi >— -— >— 



■5 2 



a 



9 o 



u 
z 



o S 

^ X 

-a 

^ "? 



Q. >, 
D X! 
O O 



o 

, 51 S 



c 
O 
U 

(N 



01 



o 



Xi 

1/1 o 

3 q 

2 " o 

is 2; q. 

C "5 "to 

2 ^ £ 

'l^ U Z! 

5; ^ -2. 

ai ^ 

i's ■'^ 

-J O -o 



PLOS ONE I www.plosone.org 



6 



May 2014 | Volume 9 | Issue 5 | e97417 



Replication Study of Osteoarthritis Candidate Genes 




PLOS ONE I www.plosone.org 



7 



May 2014 | Volume 9 | Issue 5 | e97417 



Replication Study of Osteoarthritis Candidate Genes 



< 
o 

Q 



E 
o 

Q. 

E 
>> 



VD 

■-^ rvj (N 



<— rsi 



VO LT) r— 



^ vo m o> o^ 

CO vO vO fvi m 



M3 



I 



> 
a 



tN <— -— t— 



fN rsl fN tN 



m rsl 



< 



rsi >— -— >— 



■5 2 



a 



to 
u 
z 



it 



c 

o 
u 

n 



01 



o 



S a; 



QJ ^ ^ 

C 4_i ^ 

O 3 
c 

-H o 

cn uj P 

2 a 



;1 

: CL "3 ^ 



PLOS ONE I www.plosone.org 



8 



May 2014 | Volume 9 | Issue 5 | e97417 



Replication Study of Osteoarthritis Candidate Genes 



which is particularly prevalent among women, especially post 
menopause. The rs 1799945 variant in the i/Fi? (hemochromatosis 
Fe) gene is the most common genetic factor in hereditary 
hemochromatosis. A high risk of developing OA has been 
observed in iron overload patients, who carry an HFE mutation 
[46] . Increasing evidence shows that a link between defective iron 
metabolism and tissue responses may drive the OA phenotype 
[47]. A significant increase occurs in women's iron levels during 
menopausal transition [48]. Therefore, iron increase could be a 
risk factor in age-related OA. Candidate gene association studies 
consistently showed associations between mutations within the 
HFE gene and different hand OA phenotypes [49]. Inverse 
changes have been observed between iron and estrogen levels in 
healthy women during menopausal transition [50]. 

The estrogen receptor 1 alpha (ESRl) gene has been 
investigated in several genetic studies of knee OA, but very little 
in relation to hand OA. Valdes et al. [51] reported an association 
between the gene and clinical knee OA in women but not in men. 
Bergink et al. [52] found an association between the ESRl 
haplotype and radiographic knee OA in elderly men and women. 
One study of Japanese women that investigated the association of 
PvuII (rs2234693) and Xbal (rs9340799) SNPs in the ESRl gene 
with generalized hand OA observed an approximately two-fold 
increase in the risk of severe generalized OA among participants 
carrying heterozygous genotypes for both rs2234693 and 
rs9340799 loci [53]. However, these associations have not been 
replicated in European descent for hand OA [26]; in European 
descent the ftu// variant allele has been associated with a reduced 
risk of hip OA in elderly women [54] . 

The SNPs selected for exploratory analysis in our study are 
located in potential candidate genes for OA identified by GWAS. 
We found that the variants in the BCAP29, DI02, DUS4L, DVWA, 
HLA, PARD3B, PTGS2, TGFBl and TRIBl genes were associated 
with OA across multiple sites [27-33]. We also discovered that the 
TGEBl (rsl800470, formerly known as rsl982073) variant allele 
(LeulO) was associated with a substantially increased risk of 
symptomatic DIP OA, but not radiographic OA. Transforming 
growth factor- fi (TGFB) is a multifunctional growth factor with 
widespread effects on multiple tissues, including an important role 
in cartilage matrix metabolism. Elevated TGFB 1 levels have been 
identified in the synovial fluid of patients with OA [55], and levels 
after meniscectomy are reflective of future bone and cartilage 
changes [56]. The functional role of the LeulOPro polymorphism 
remains largely unknown, but some evidence exists that it affects 
TGFBl secretion and functions in hepatic cells [57]. Earlier, this 
polymorphism was believed to be associated with a radiographic 
spinal OA among Japanese women [32]. 

Our study had both strengths and limitations. The main 
limitation was the relatively small sample size, due to which the 
study had low power to detect small effects. One major strength 
was that it consisted of random samples of ethnically relatively 
homogenous Finnish origin. The Finnish population is one of the 
best-studied genetic isolates, which originated from a small 
founder population some 2000 years ago. Therefore, the Finnish 
population has a relatively homogenous gene pool [58] offering 
optimal material for association studies. We also chose to have 
only women in our study, thus findings that might be restricted to 
men were left undetected. Moreover, the age range (45-63 years) 
of the participants was selected to cover the age of occupationaUy 
active women, whose prevalence of hand OA is rapidly increasing. 
Finally, dentists and teachers have a similar level of education, 
although occupational exposures related to hand use, and lifestyle 
factors such as obesity, differed. 
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The current findings were also unlikely to have been influenced 
by selection bias. The prevalence of the studied hand OA 
phenotypes was similar to those observed in other studies [59,60]. 
Yet, the 'healthy worker efiect' may have led to an underestima- 
tion of associations between SNPs and hand OA, especially for 
symptomatic DIP OA. Dentists suffering from severe pain in the 
hand joints may not be able to continue in their profession and 
therefore were outside the occupationaUy active group that was 
our target population. In contrast, teachers with symptomatic 
hand OA may remain active in their profession. Indeed, 
occupation-stratified analyses revealed that most of the studied 
SNPs had a different efiect on hand OA in the two occupational 
groups. 

To summarize, by replicating and pardy confirming earlier 
studies of OA susceptibility polymorphisms, our results further 

support the theory that the A2BP1 and TBGFl genes are hand OA 
susceptibility genes. OA is a globally significant medical, social, 
and economic problem and therefore merits attention in order to 
develop better prevention methods and therapies that can be 
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